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By combining data from seven cosmic microwave background experiments (including the latest 
WMAP results) with large scale structure data, the Hubble parameter measurement from the Hubble 
Space Telescope and luminosity measurements of Type la supernovae we demonstrate the bounds 
on the dark energy equation of state wq to be —1.38 < wq < —0.82 at the 95% confidence level. 
Although our limit on wq is improved with respect to previous analyses, cosmological data does not 
rule out the possibility that the equation of state parameter wq of the dark energy Q is less than -1. 
We present a tracking model that ensures wq < — 1 at recent times and discuss the observational 
consequences. 


I. INTRODUCTION. 

There is now a growing body of evidence that the 
evolution of the universe may be dominated by a dark 
energy component Q^ith present-day energy density 
fraction UIq ~ 2/3 jlj. Although a true cosmologi¬ 
cal constant A may be responsible for the data, it is 
also possible that a dynamical mechanism is at work. 
One candidate to explain the observations is a slowly- 
rolling dynamical scalar “quintessence” field |2|-0. 
Another possibility, known as “k-essence” 011 , is 
a scalar field with non-canonical kinetic terms in the 
lagrangian. Dynamical dark energy models such as 
these, and others @-Q have an equation of state, 
wq = Pq/pq which varies with time compared to 
that of a cosmological constant which remains fixed 
at wq- a = — 1. Thus, observationally distinguishing 
a time variation in the equation of state or finding 
wq different from —1 will rule out a pure cosmolog¬ 
ical constant as an explanation for the data, but be 
consistent with a dynamical solution. 

In the past years many analyses of several cosmolog¬ 
ical datasets have been produced in order to constrain 
wq (see e.g. E2 and references therein). In these 
analyses the case of a constant-with-redshift wq, in 
the range wq > — 1 was considered. The assumption 
of a constant wq is based on several considerations: 
first of all, since both the luminosities and angular 
distances (that are the fundamental cosmological ob¬ 
servables) depend on wq through multiple integrals, 
they are not particularly sensitive to variations of wq 
with redshift (see e.g. [Dj, (3)- Therefore, with cur¬ 
rent data, no strong constraints can be placed on the 
redshift-dependence of wq. Second, for most of the 
dynamical models on the market, the assumption of a 
piecewise-constant equation of state is a good approx¬ 
imation for an unbiased determination of the effective 
equation of state fll3|) 

f w Q (a)fl Q (a)da 

*eff ~ / w ( 1 ) 

J \lQ(a)aa 

predicted by the model. Hence, if the present data is 
compatible with a constant wq = — 1 , it may not be 


possible to discriminate between a cosmological con¬ 
stant and a dynamical dark energy model. 

The limitation to wq > —1, on the contrary, is a 
theoretical consideration motivated, for example, by 
imposing on matter (for positive energy densities) the 
null energy condition, which states that T^N^N 1 * > 
0 for all null 4-vectors iV M . Such energy conditions 
are often demanded in order to ensure stability of the 
theory. However, theoretical attempts to obtain wq < 

— 1 have been considered 0 EH EH EH EH while a 
careful analysis of their potential instabilities has been 
performed in El- 

Moreover, Maor et al. {3 have recently shown that 
one may construct a model with a specific z-dependent 
wq{z) > —1, in which the assumption of constant wq 
in the analysis can lead to an estimated value w e s < 

— 1. This further illustrates the necessity of extending 
dark energy analyses to values of wq < — 1. 

In this paper we combine constraints from a variety 
of observational data to determine the currently al¬ 
lowed range of values for the dark energy equation of 
state parameter wq. The data used here comes from 
six recent Cosmic Microwave Background (CMB) ex¬ 
periments, from the power spectrum of large scale 
structure (LSS) in the 2dF 100k galaxy redshift sur¬ 
vey, from luminosity measurements of Type la super¬ 
novae (SN-Ia) El and from the Hubble Space Tele¬ 
scope (HST) measurements of the Hubble parameter. 

Our analysis method and our datasets are very sim¬ 
ilar to the one used in a recent work by Hannestad 
and Mortsell (01)- We will compare our results with 
those derived in this earlier paper in the conclusions. 

In the next section, we demonstrate the plausibil¬ 
ity of wq < — 1 by presenting a class of theoretical 
models in which this result may be obtained explic¬ 
itly. In our model the equation of state parameter is 
approximately piecewise constant and hence provides 
a specific example of a model which would fit the data 
described in the remainder of the paper. The meth¬ 
ods used to obtain combined constraints on the dark 
energy equation of state are described in section mu 
Our likelihood analysis is presented in section ITvl and 
our summary and conclusions are given in section Ivl 
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II. A MODEL WITH wq < 1. 

It is a simple exercise to show that a conventional 
scalar field lagrangian density cannot yield an equa¬ 
tion of state parameter wq < — 1 . There are, how¬ 
ever, a number of ways in which the lagrangian can 
be modified to make wq < — 1 possible. For example, 
one may reverse the sign of the kinetic terms, leading 
to interesting cosmological and particle physics behav¬ 
ior Ei-ta. 

Let us motivate the study of wq < — 1 cosmologies 
by describing a class of models, with positive kinetic 
energy , in which such evolution arises. The model we 
consider is very much in the spirit of k-essence jUfT]- 
Comments on the similarities and differences between 
the two are briefly discussed at the end of this section. 

Consider a theory of a real scalar field </>, assumed to 
be homogeneous, with a non-canonical kinetic energy 
term. The lagrangian density is 

C = mg{X) - V(</>) , (2) 

where g(X ) are positive semi-definite functions, 
V(4>) is a potential and X = </> 2 /2. The energy- 
momentum tensor for this field is straightforward to 
calculate and yields the usual perfect fluid form with 
pressure p and energy density p given by 

V = C = mg(X) - V(0) , (3) 


P = 


2X MX} {x) 

dX ; 


m+vw). ( 4 ) 


Thus, defining w$ = pip one obtains 

a{x)m - v{</>) 


W<t> = 


2 X&p--g{X) m+Vtf) 


(5) 


If the lagrangian 0 is to yield < — 1 then 0 
implies 


for a detailed stability analysis of models with wq < 
— 1). This yields 


2 = dp = ppi_ = g'(X) 

° s ~ dp p.x g'{X) + 2 Xg»(X) 


where a subscript , X denotes a partial derivative with 
respect to X. Since we have already specified g'(X) < 
0 this may be written as 


ff"(X) < 


g'(x) 

2X 


(9) 


Notice the difference between this class of models 
and the k-essence family in terms of the potential V (</>) 
and the constraints placed on the functions g(X) and 
/(</>)■ 

Let us illustrate these constraints with a simple ex¬ 
ample g{X) = e ~ aX , with a > 0. This function 
trivially satisfies the constraints The con¬ 

straint 0 then yields 


m 


> (2aX + l)e~ aX . 


( 10 ) 


In the asymptotic regions this becomes V/f > 0 as 
X —> oo and V/ f > lasX—>0. This may be satisfied 
without a particularly large potential by arranging an 
appropriately behaved /(</>), since 0 constrains only 
the ratio V/f, making fine-tuning issues less severe. 

Let us now assume a (flat) Friedmann, Robertson- 
Walker (FRW) ansatz for the space-time metric 


ds 2 = — dt 2 + a{t) 2 dHt 2 , (11) 


with a(t) the scale factor. The resulting Einstein equa¬ 
tions then become the Friedmann equation 



and the acceleration equation 


( 12 ) 


g\X) < 0 , (6) 

where g’{X) = dg(X)/dX and we have used > 0 
and X > 0. We therefore require that g(X) be a 
strictly monotonically decreasing function. It is inter¬ 
esting to note in passing that, provided is positive 
semi-definite, the functional forms of and V((f>) 
play no role in determining whether w</, is less than or 
greater than -1. 

However, a constraint involving the potential does 
arise from the requirement that the energy density of 
the theory satisfy p > 0. This yields 

g(X) - 2Xg\X) <J^- (7) 

A necessary condition that the theory be stable is 
that the speed of sound of </> be positive [U (see 0 


a 47 tG . „ . 

- =-5“(p + 3p). (13) 

a 3 

One then solves these equations along with those for 
the scalar field. In the case of k-essence with wq > —1 
it has been shown 0-0 that tracking behavior can 
be obtained. This means that for a wide range of 
initial conditions, the energy density of the field <fi 
naturally evolves so as to track the energy density in 
matter, providing some insight into why dark energy 
domination began only recently in cosmic history. In 
our model, in which we have included a potential for 
(f> and are in the regime wq < — 1, the analysis be¬ 
comes somewhat more involved. Nevertheless, it can 
be shown that tracking behavior persists. However, as 
in all rolling scalar models, some fine-tuning remains, 
since one must ensure the right amount of dark energy 
density today. 















III. COMPARISON WITH OBSERVATIONS: 
METHOD 
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We restrict our analysis to flat models, for which 
the effects of dark energy with wq > -1 on the an¬ 
gular power spectrum of the CMB anisotropies have 
been carefully analyzed (see e.g. JJj and references 
therein). The main effect of changing the value of wq 
on the CMB anisotropies is to introduce a shift by a 
linear factor TZ in the Cspace positions of the acoustic 
peaks in the angular power spectrum p5 |. This shift 
is given by 


n = yf(l^n Q )y , (14) 


where 

[ Zdec dz 

V = / , ■ 

Jo ^/( 1 - Oq )(1 + z ) 3 + O q (1 + 2 ) 3 ( 1 +- o ) 

(15) 

In order to illustrate this effect, we plot in Figure^a 
set of theoretical power spectra, computed assuming a 
standard cosmological model with the relative density 
in cold dark matter fl c d m h 2 = 0.12, that in baryons 
n b h 2 = 0.022, with Hubble parameter h = 0.69 but 
with wq varied in the range (—4,—0.5). It is clear 
that decreasing wq shifts the power spectrum towards 
smaller angular scales 6 ~ l _1 . 

In considering the CMB power spectrum, it is im¬ 
portant to note that there is some degeneracy among 
the possible choices of cosmological parameters. First 
of all, the shift produced by a change in wq can easily 
be compensated by a change in the curvature. How¬ 
ever, degeneracies still exist even when we restrict to 
flat models. To emphasize this we plot in Figure [3 
some degenerate spectra, obtained by keeping Q b h 2 , 
VlMh 2 , and TZ fixed, in a flat universe. In practice, 
in order to preserve the shape of the spectrum while 
decreasing wq , one has to increase Qq . For flat mod¬ 
els, one must therefore decrease VLm and, since Uifli 2 
must be constant, increase h. Therefore, even if the 
CMB spectra are degenerate, combining the CMB in¬ 
formation with priors on Am and h can be extremely 
helpful in bounding wq. 

On large angular scales the time-varying New¬ 
tonian potential after decoupling generates CMB 
anisotropies through the Integrated Sachs-Wolfe 
(ISW) effect. This is clearly seen in Figure |3 and is 
more pronounced for more negative values of wq . The 
effect depends not only on the value of wq but also on 
its variation with redshift. However, this is difficult to 
disentangle from other cosmological effects. 

In all our analysis we will neglect perturbations in 
the dark energy component. The reasons for this sim¬ 
plification are twofold: on the one hand we prefer our 
analysis to remain as model-independent as possible, 
so that the results obtained here are not affected by 
the choice of a particular dark energy model. The 



FIG. 1: The effect of varying wq on the COBE-normalized 
CMB angular power spectrum and present CMB data. 
Since the shift of the power spectrum is proportional to 
7 Z£, the discrepancy is more important for higher values of 
l. 



FIG. 2: Degenerate CMB power spectra The models are 
computed assuming flatness (Qtot = 9m + Hq = 1). 
On large angular scales the Integrated Sachs Wolfe effect 
breaks the degeneracy for highly negative values of wq . In 
general, the degeneracy of the spectra can be broken with 
a strong prior on h or on 9 m. 


study of the perturbations in particular quintessen¬ 
tial models goes beyond the scope of this paper. On 
the other hand, this approximation is also satisfied 
in a broad class of models and it is not completely 
straightforward to conclude that the inclusions of per¬ 
turbations, while consistent with General Relativity, 
would yield a better approximation for a particular 
model of dark energy. As an example, in Figure^ we 
plot the CMB power spectra for w = —0.8 computed 
with and without assuming adiabatic perturbations 
in the dark energy fluid as in (2^|. As we can see, 
since dark energy dominates the overall density only 
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FIG. 3: Effect of including perturbations in the dark en¬ 
ergy fluid with a constant equation of state in the CMB 
power spectrum. Since dark energy dominates at late red- 
shifts, the effect is present only on large angular scales. 


well after recombination, the major effects are only 
on very large scales. We found that the inclusion of 
perturbations has no relevant effect on our results for 
models with w > — 1 where the computations of the 
perturbations are meaningful. 

In order to bound wq, we consider a template of 
flat, adiabatic, Q-CDM models computed with CMB- 
FAST H. We sample the relevant parameters as fol¬ 
lows: f Icdmh 2 = 0.01, ...0.40, in steps of 0.01; f Ibh 2 = 
0.001,..., 0.040, in steps of 0.001, Hq = 0.0,..., 0.95, in 
steps of 0.05 and wq = —3.0,..., —0.4 in steps of 0.05. 
Note that, once we have fixed these parameters, the 
value of the Hubble constant is not an independent 
parameter, since it is determined through the flatness 
condition. We adopt the conservative top-hat bound 
0.45 <h< 0.85. 

We allow for a reionization of the intergalactic 
medium by varying the Compton optical depth pa¬ 
rameter t c over the range r c = 0.05,..., 0.40 in steps 
of 0.05. 

For the CMB data we use the recent tempera¬ 
ture and cross polarization results from the WMAP 
satellite ([Hj) using the method explained in (f2t|) 
and the publicly available code on the LAMBDA 
web site. We further include the results from the 
BOOMERanG-98 [p. DASI [H, MAXIMA-1 H, 
CBI [3Cj, VSAE l3l|. and Archeops [111 experiments 
by using the publicly available correlation matrices 
and window functions. We consider 7%, 10%, 4%, 
5%, 3.5% and 5% Gaussian distributed calibration 
errors for the Archeops, BOOMERanG-98, DASI, 
MAXIMA-1, VSA, and CBI experiments respectively 
and include the beam uncertainties using the analyt¬ 
ical marginalization method presented in [33ll . The 
likelihood £ for a given theoretical model is defined 

by 

- 2 In £ = (Cg - C%)M bb , (C# - C B ,) , (16) 

where M BB i is the Gaussian curvature of the like¬ 


lihood matrix at the peak and C B is the theore- 
ical/experimental signal in the bin (|34j). 

In addition to the CMB data we also consider the 
real-space power spectrum of galaxies in the 2dF 100k 
galaxy redshift survey using the data and window 
functions of the analysis of Tegmark et al. (35]. To 
compute the likelihood function C 2dF for the 2dF sur¬ 
vey we evaluate pi = P{ki), where P(k) is the theoret¬ 
ical matter power spectrum and ki are the 49 k-values 
of the measurements in |3jj ■ Therefore 


— 2 In C 2dF 


E 


[Pj - ( Wp ),;] 2 

dP 2 


(17) 


where Pi and dPi are the measurements and corre¬ 
sponding error bars and W is the reported 27 x 49 
window matrix. We restrict the analysis to a range of 
scales over which the fluctuations are assumed to be 
in the linear regime (k < 0.1/i -1 Mpc). When combin¬ 
ing with the CMB data, we marginalize over a bias b 
considered to be an additional free parameter. 

We also incorporate constraints obtained from the 
luminosity measurements of Type la supernovae (SN- 
Ia). In doing this note that the observed appar¬ 
ent bolometric luminosity m B is related to the lu¬ 
minosity distance g(l, measured in Mpc, by m B = 
M + 5 logdi(z) + 25, where M is the absolute bolomet¬ 
ric magnitude. Note also that the luminosity distance 
is sensitive to the cosmological evolution through an 
integral dependence on the Hubble factor 


d L = (l + z) [ 

Jo 


dz' 


H(z',SlQ,Q, M ,w q ) 


(18) 


We evaluate the likelihoods assuming a constant equa¬ 
tion of state, such that 


H 2 (z) = Hq J2 fia(l + *) (3+3wo) , (19) 

a 


where the subscript a labels different components of 
the cosmological energy budget. The luminosity m e s 
predicted from the observations is then calculated 
by calibration with low-z supernovae observations for 
which the Hubble relation d^ « Hqcz is obeyed. We 
calculate the likelihood £ SN using the relation 


£ sn 


£o exp 


X 2 (H qMm,wq) 
2 


( 20 ) 


where Cq is an arbitrary normalisation and y 2 is eval¬ 
uated using the observations of ijj and marginalising 
over H 0 . Finally, we also consider the lcr contraint on 
the Hubble parameter, h = 0.71 ±0.07, obtained from 
Hubble Space Telescope (HST) measurements [36]. 


IV. COMPARISON WITH OBSERVATIONS: 
RESULTS 

Table I shows the 2-cr constraints on wq in a flat 
universe for different combinations of priors, obtained 
after marginalizing over all remaining parameters. 
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CMB+HST -1.65 < wq < -0.54 

0.19 < Om < 0.43 
CMB+HST+BBN -1.61 < wq < -0.57 

0.20 < S7 M < 0.42 
CMB+HST+SN-Ia -1.45 < wq < -0.74 

0.21 < Om < 0.36 
CMB+HST+SN-Ia+2dF -1.38 < w Q < -0.82 
0.22 < fi M < 0.35 


TABLE I: Constraints on wq and f Im = 1 — Hq using 
different priors and datasets. We always assume flatness 
and that the age of the universe to > 10 Gyr. The 20- 
limits are found from the 2.5% and 97.5% integrals of the 
marginalized likelihood. The HST prior is h = 0.71 ±0.07, 
while for the BBN prior we use the conservative bound 
Q b h 2 = 0.020 ± 0.005. 


It is clear that wq is poorly constrained from CMB 
data alone, even when the strong prior on the Hubble 
parameter from HST, h = 0.71 ± 0.07, is assumed. 
Adding a Big Bang Nucleosynthesis prior, fl b h 2 = 
0.020 ± 0.005, has a small effect on the CMB±HST 
result. Adding SN-Ia data breaks the CMB — wq 
degeneracy and improves the limits on wq, yielding 

— 1.45 < wq < —0.74. Finally, including data from 
the 2dF survey further breaks the degeneracy, giving 

— 1.38 < wq < —0.82 at 2-a. Also reported in Ta¬ 
ble I are the constraints on CIm- The combined data 
suggests the presence of dark energy with high signif¬ 
icance, even if one only considers CMB±HST data. 

It is interesting to project our likelihood onto the 
(Om, wq) plane. In Figure U| we plot the likelihood 
contours in the (Om, wq) plane from our joint analy¬ 
ses of CMB±SN-Ia±HST±2dF data. As we can see, 
there is strong supporting evidence for dark energy. 
A cosmological constant with wq = — 1 is in good 
agreement with all the data and the most recent CMB 
results improve the constraints from previous and sim¬ 
ilar analyses (see e.g. |2Cjh 


V. CONCLUSIONS 

In this paper we have provided new constraints 
on the dark energy equation of state parameter wq 
by combining recent cosmological datasets. We find 
— 1.38 < wq < —0.82 at the 95% confidence level, 
with best-fit model Hm = 0.27 and wq = —1.05. A 
cosmological constant is a good fit to the data. When 
comparison is possible (i.e. restricting to similar priors 
and datasets), our analysis is compatible with other 
recent analyses of wq (e.g. see a 0 and refer¬ 
ences therein), however, our lower bound on wq is 
much tighter than the one recently reported in l20ll . In 
particular, we found that the CMB±HST dataset can 
already provide an interesting lower limit on wq , while 



n matter 


FIG. 4: Likelihood contours in the (I2 m, wq) plane for 
the joint CMB±HST±SN-Ia±2dF analysis described in 
the text. We take the best-fit values for the remaining pa¬ 
rameters. The contours correspond to 0.32, 0.05 and 0.01 
of the peak value of the likelihood, which are the 68%, 
95% and 99% confidence levels respectively. Also plotted 
are the likelihood contours from type-la Supernovae alone. 

in f2Qj no constraint was obtained. Part of the discrep¬ 
ancy can be explained by our updated CMB dataset 
with the new Archeops, Boomerang, CBI, VSA and, 
mostly, WMAP results. However, our CMB power 
spectra in Fig.l are in disagreement with the same 
spectra plotted in the Fig.2 of J2(j where the depen¬ 
dence on wq seems limited only to the large-scale ISW 
term. 

In the range w > — 1 our results are in very good 
agreement with those reported by Spergel et al. l37l . 
which is using a different analysis method based on 
a Monte Carlo Markov Chain and a slightly different 
CMB dataset. 

We found that including models with w < — 1 does 
not significantly affect the results obtained under the 
assumption of w > — 1. In this respect, our findings 
are a useful complement to those presented in [37]. 

As in P and in most of previous similar analysis 
the constraints obtained here have been obtained un¬ 
der several assumptions: the equation of state is red- 
shift independent, the perturbations in the dark en¬ 
ergy fluid are negligible and/or its sound speed never 
differs from unity in a significant way. It is impor¬ 
tant to note that our result apply only to models well 
described by these approximations. 

We have also demonstrated that, even by apply¬ 
ing the most current constraints on the dark energy 
equation of state parameter wq , there is much uncer¬ 
tainty in its value. Interestingly, there is a distinct 
possibility that it may lie in the theoretically under¬ 
explored region wq < — 1. To illustrate this we have 
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provided a specific model in which wq < — 1 is at¬ 
tained, and which satisfies the assumption that wq is 
approximately piecewise constant, as used in the data 
analysis. An observation of a component to the cosmic 
energy budget with iuq < — 1 would naturally have 
significant implications for fundamental physics. Fur¬ 
ther, depending on the asymptotic evolution of wq, 
the fate of the observable universe MM may be 
dramatically altered, perhaps resulting in an instabil¬ 
ity of the spacetime jig] or a future singularity. 

If we are to understand definitively whether dark 
energy is dynamical, and if so, whether it is consis¬ 
tent with wq less than or greater than —1, we will 
need to bring the full array of cosmological techniques 
to bear on the problem. An important contribu¬ 
tion to this effort will be provided by direct searches 
for supernovae at both intermediate and high red- 
shifts |43|. Other, ground-based observations (43 will 
allow complementary analyses, including weak grav¬ 
itational lensing [44| and large scale structure sur¬ 
veys HH to be performed. 

At present, however, while the data remain con¬ 


sistent with both a pure cosmological constant A, and 
with dynamic classes of models mm , nature may 
be telling us that the universe is an even stranger place 
than we had imagined. 
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